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Abstract. Heteronuclear alkali-metal dimers represent the class of molecules of 
choice for creating samples of ultracold molecules exhibiting an intrinsic large 
permanent electric dipole moment. Among them, the KCs molecule, with a permanent 
dipole moment of 1.92 Debye still remains to be observed in ultracold conditions. Based 
on spectroscopic studies available in the literature completed by accurate quantum 
chemistry calculations, we propose several optical coherent schemes to create ultracold 
bosonic and fermionic KCs molecules in their absolute rovibrational ground level, 
starting from a weakly bound level of their electronic ground state manifold. The 
processes rely on the existence of convenient electronically excited states allowing an 
efficient stimulated Raman adiabatic transfer of the level population. 


Prospects for the formation of ultracold polar ground state KCs molecules via an optical process2 


PACS numbers: 33.20.-t,37.10.Mn,37.10.Vz,34.20.-b 



Prospects for the formation of ultracold polar ground state KCs molecules via an optical process?) 

1. Introduction 

Dilute atomic and molecular gases at ultracold temperatures (T = E/ks 
1 millikelvin) offer the fascinating opportunity of long observation time allowing for 
measurements with unprecedented accuracy. For instance, due to their extremely low 
relative velocity in such gases, particles have their maximal presence probability at 
large mutual distances R, well beyond the range of electron exchange. Therefore the 
dynamics of ultracold gases is dominated by their weak long-range (van der Waals) 
interaction varying as which is isotropic for identical particles in free space with 
spherical symmetry (like atoms, or molecules with vanishing total angular momentum). 

Once the particles are immersed in an external magnetic or electric field, their intrinsic 
properties (permanent magnetic or electric dipole moment) induce the anisotropy of 
their long-range interaction, which can now vary as R~^ and depends on the relative 
orientation of their molecular axis, or of their angular momentum [1, 2]. Manifestations 
of anisotropic interactions have already been observed experimentally with ultracold 
quantum degenerate gases of magnetic atoms [3, 4, 5, 6], and during ultracold collisions 
between KRb polar molecules (he. possessing a permanent electric dipole moment 
in their own frame) [7, 8]. Such so-called ultracold dipolar gases are expected to 
reveal novel physical phenomena for instance in the context of quantum degeneracy 
where hamiltonians involved in condensed matter physics could be simulated with the 
opportunity for controlling the interaction between particles with external fields (see 
for instance the review papers of Refs. [9, 10]). As for molecules, the recent review 
articles of Refs. [11, 12] provide in-depth presentations of theory and experiments of 
collisions and reactions with ultracold molecules, emphasizing on their implications in 
the development of the new research area of ultracold chemistry dominated by quantum 
mechanical effects [13]. 

Despite amazing experimental progress, the main challenge for experimentalists 
dealing with ultracold molecules is still their formation as a gaseous sample with 
sufficient number density and with a good control of their internal state. An overview 
of the various methodologies to create ultracold neutral molecules and of their potential 
opportunities and applications is available in several review articles [14, 15, 16], and we 
will not cover them here. In brief, there are two classes of approaches to obtain ultracold 
ground state molecules: (i) manipulating pre-existing polar diatomic or polyatomic 
molecules with external magnetic or electric fields to design slow molecular beams 
[17, 18, 19], and for some specific polar species, cooling diatomic molecules with laser 
[20, 21, 22, 23, 24]; (ii) associating a pair of ultracold atoms into an ultracold molecule 
using laser photoassociation (PA) toward an electronic excited state followed by radiative 
emission (RE) to the ground state [25, 26], or magnetoasssociation (MA) [27] in a 
weakly bound level of the electronic ground state manifold via magnetically tunable 
Feshbach resonances [28], with a subsequent stimulated radiative transfer (SRT) process 
to populate the lowest molecular bound level of the ground state. 

We focus for the rest of this paper on the latter option (MA-I-SRT), which has 
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been successfully demonstrated experimentally in a still limited number of cases with 
homonuclear molecules like Cs 2 [29] and Rb 2 [30], and heteronuclear molecules like KRb 
[31] and RbCs [32, 33]. This is a quite general method for the class of alkali-metal 
diatomic molecules which all possess a wealth of Feshbach resonances (see Ref. [34] for 
KCs isotopologues). The population transfer toward the lowest bound level of the ground 
state is achieved via the well-known coherent process of stimulated Raman Adiabatic 
Passage (STIRAP, [35, 36, 37]). The efficiency of the transfer relies on the identihcation 
of a pair of so-called pump and dump electric dipole allowed transitions with comparable 
Rabi frequencies, and thus on the detailed knowledge of the spectroscopy of the molecule 
of interest. In a previous paper [38], hereafter referred to as paper I, we modeled the 
STIRAP approach for the bosonic and fermionic KRb molecules. Using up-to-date 
spectroscopic data, we analyzed the efficiency of several transition schemes over the 
entire range of accessible laser frequencies determined by the excited electronic states, 
and we confirmed the suitability of the experimentally chosen scheme for KRb. We also 
investigated STIRAP efficiency on a limited range of frequencies for RbCs [39, 32]. 

Our study concerns the formation of ultracold bosonic ^^K^^^Cs and fermionic 
40p^i33Qg polar molecules in their absolute ground state by STIRAP, which has not 
yet been achieved experimentally. The present investigation aims at guiding ongoing 
experiments in the choice of their laser set-up to implement the STIRAP scheme. Its 
main outcome is that the most efficient STIRAP schemes do not generally correspond to 
the intuitive picture delivered by the Franck-Condon principle describing the strongest 
molecular transitions in terms limited to the spatial overlap of vibrational wave 
functions. The KCs species exhibits several specific properties which contrast with 
the KRb one studied in paper I. It possesses an intrinsic permanent electric dipole 
moment (PEDM) of 1.92 D almost 4 times larger than the KRb one [40]. The potential 
well depth of the KCs molecular ground state being larger than the Cs 2 one, but 
smaller than the K 2 one, the KCs molecule is stable in its lowest ground state level 
against ultracold collisions with surrounding Cs atoms and other KCs molecules [41], 
but not with surrounding K atoms. These features will represent a decisive advantage 
for their further manipulation to create for instance a quantum degenerate molecular 
gas without the request of trapping molecules inside an optical lattice. Moreover, the 
KCs spectroscopy has been already quite well investigated experimentally for the 
ground state and the lowest triplet state a^S+ [42, 43], and for several excited states 
[44, 45, 46, 47, 48, 49]. The spectrum of KCs Feshbach resonances based on these results 
has also been modeled in detail [34]. 

The present paper is organized as follows. We first recall in Section 2 the basic 
principle of STIRAP, and we characterize the initial, and final molecular electronic 
states chosen for the implementation of STIRAP in ^®KCs and in '^^KCs (the index for 
the Cs mass will be omitted in the rest of the paper). The choice of the intermediate 
state for STIRAP is discussed in Section 3, emphasizing on the necessary knowledge of 
the requested molecular structure data, i.e. potential energy curves (PECs), transition 
electric dipole moments (TEDMs) and spin-orbit couplings (SOCs). Three transition 
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paths for STIRAP in KCs relying on different kinds of couplings between molecular 
levels are identified and their efficiencies are compared among each other (Section 4). 
Experimental prospects are discussed in Section 5 in the perspective of the extension of 
such studies to other alkali-metal polar diatomic species. 

When appropriate, the atomic unit of length (1 oq = 0.052917721092 nm) and of 
dipole moment (1 a.u.= coq = 2.541 580 59 D) will be used. 

2. Model for STIRAP with KCs molecules 

The STIRAP principle has been proposed by Bergmann et a/. [35, 36] and further 
discussed for instance in Ref. [37] in the context of coherent photoassociation of ultracold 
atoms. We also presented a summary in paper I and we only recall here a few aspects 
which are relevant for the present study. 

The central idea of STIRAP is to adiabatically transfer the population of a quantum 
system from an initial state | i) to a well-dehned hnal state | g) via an intermediate 
excited state | e), in such a way that the | e) state is actually not populated. This 
is achieved by cleverly shaping two laser pulses overlapping in time inducing the | i) 

^ I e) (pump) transition and the | e) \ g) (dump) transition. The Hamiltonian of 

the system dressed by the pulses admits a ’’‘dark’” eigenstate which cannot radiatively 
decay and insures the transfer from | i) to | ) without loss of population. One can show 
that the optimal efficiency of the transfer is reached when one hnds a level | e) such as 
the amplitude of the time-dependent Rabi frequencies 

Oei ( e I d ■ T/punjp 10/^1 { 9 \ ^ ' A'dump I ^ ^ (f) 

for the pump and dump transitions are equal. In Eq.(l) d{R) is the electronic transition 
dipole moment function (TEDM), and Epu^p (resp. Ejump) is the amplitude of the 
laser held driving the pump (resp. dump) transition, associated to an intensity /pump 
(resp /dump)- This equality is achieved either with equal TEDM matrix elements for the 
pump and dump transitions assuming equal intensities for both pulses which is often 
convenient in experimental setups, or by slight adjustments of the laser intensities within 
the experimental feasibility to approach the strict equality of the Rabi frequencies. 

Here | ( 7 ) is the lowest level ux = 0 of the electronic ground state of KCs, and 

I i) is a weakly-bound level of the and a^E+ state manifolds (hereafter referred 

to as the X and a states) coupled by the hyperhne interaction. Figures 1 and 2 display 
two possible choices of excited electronic state for the intermediate state | e) that are 
discussed below. 

The weakly-bound initial state | i) results from the magnetoassociation of a pair of 
ultracold ^®K (or “^^K) and Cs atoms by tuning an external magnetic held onto a Feshbach 
resonance of the pair. Such a so-called Feshbach molecule is populated in a high-lying 
rovibrational level with a combination of triplet {S = 1) and singlet {S = 0) characters. 

The mixing coefficients depend on the choice of the Feshbach resonance, experimentally 
investigated in Ref. [43] and accurately modeled in Ref. [34]. For a molecular system. 
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there is an additional reqnirement for STIRAP to work, that the radial wavefnnctions 
of the I i) and | e) on one hand, and of the | e) and | g) levels on the other hand should 
overlap each other in a region where the TEDMs are not vanishingly small. As it can 
be seen from Figs, la and 2a, this is well achieved if the | i ) levels contains a significant 
component on the a^S’*' state which has an inner classical turning point in the suitable 
range of R. Just like in paper I, the first hypothesis of our model is to choose | i) with a 
pure triplet character, so that one can describe it with a radial wavefunction belonging 
to the single a^S+ PEC with approximately the same binding energy as the one of the 
Feshbach molecule bound level. This is the case of the uppermost a^S’*' level assigned to 
Va = 35 with a binding energy of about 0.05 GHz at zero magnetic field, which possesses 
a triplet character up to 95% [43]. Note, however, that a pump transition starting from 
a pure singlet level would, in principle, be possible (see the dashed arrows in Fig. la) 
but is not discussed further here. 

Such an assumption is reasonable, as even if there is a significant singlet component 
in the chosen level, the nodal structure of the triplet component of the radial 
wavefunction will not be affected in the region of the inner turning point of the a^E’*' 

PEC, while its amplitude may be changed. Therefore the matrix elements in Eq. (1) 
would be affected only through a global scaling factor. An example of such a wave 
function concerns a level resulting from the mixture of two a^E’*' levels {va = 32 and 
Va = 33) and two singlet levels {vx = 102 and Va = 103) (see Fig.6 of Ref.[43]). In a 
recent proposal, Klincare et al. [48] proposed a STIRAP implementation based on a 
pump transition mainly acting around the outer turning point of the X PEC, using a 
similar hypothesis of a pure singlet weakly-bound level as the | i) state. 

3. The choice of the intermediate STIRAP levels in KCs 

The |e) state must be optically coupled to both |i) =| a^E+ Va = 35) and l^?) = 

I A^E’*' ux = 0) and thus must exhibit favorable transition dipole moments and good 
spatial overlaps with | i) and | g) vibrational wavefnnctions. As in paper I, such a 
mixed singlet/triplet character is offered by the excited electronic states converging to 
the K(4s) -|- Cs(6p) dissociation limit, namely b^H, A^E"*", c^E"*" and B^H, hereafter 
referred to as b. A, c and B states respectively (see Figs. la and 2a). These states 
are significantly affected by the spin-orbit (SO) interaction, resulting in coupled states 
with both spin characters, labeled with the Hund’s case (c) quantum number C = O’*", 1 
for the projection of the total electronic angular momentum on the molecular axis. 

In order to provide predictions to the experimentalists, it is crucial to rely on all 
available spectroscopic information about these states, so that the corresponding PECs 
are built piecewise, combining spectroscopic and quantum chemistry determinations. 

Note that the implementation proposed in Ref. [48] relies on a higher electronic excited 
state, the (4)^E''' state correlated to the K(4s)-|-Cs(5(J) dissociation limit, perturbed by 
neighboring triplet states. 

Unlike the KRb case, our study of the C = O’*" symmetry has been greatly facilitated 
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by the extensive spectroscopic study from Refs. [45, 44, 49], providing the relevant 
PECs (Fig. la) and i?-dependent spin-orbit couplings (SOC) (Fig. lb). Following these 
authors, we use a four-coupled-channel model which accounts for the dominant SO 
interaction between the state and the O = O’*" component bo of the b^II state, 

as well as the interaction with other molecular states responsible for the asymmetric 
splitting between the O = 0,1,2 components of the b^II state. In addition it includes 
the weak rotational interactions with the other O = 1,2 components bi and b 2 of the 
b^n state, scaling with [3 = hf /{2p,R‘^) where /r is the KCs reduced mass. The resulting 
R-dependent potential energy matrix is expressed, for a given total angular momentum 
J (including the total electronic angular momentum and the rotation of the molecule, 
but not the nuclear spins, and X = J{J -|- 1)), as the sum of Born-Oppenheimer (BO) 

PEC matrix Vbo and SOC matrix ITsl? ^ 


Tbo + = 

I 62) 

/ V'^{R) + A+{R) 

-/3(l-76)\/2(2t-2) 
0 

V 0 


|61) 

-P{l-^b)V2{X-2) 

Ki(^) 

-/3(l-7b)^ 


I 60) 


0 

-/3(1-76)V^ 
no(^) - d-(it) 
-V 2 C°(i?) 


( 2 ) 

Id) 

-PCAblV^ 

-V 2 C°(i?) 

Vf{R) 


For compactness purpose we used the notations for the potential energies including the 
centrifugal term: Vf^{R) = H(R)+/?(!- 7 fe)(X-2), V',{R) = H(^)+/9(1-7b)(^ + 2), 
14'o(R) = 14(R)+/?(1-76)(X + 2), andVfiR) = 14(R)+ /9(1-7b)(^ + 2). This matrix 
reduces to a 3 x 3 form if J = 1. Note that, strictly speaking, the labels ”‘so”’ and O’*" 
in hPsl? \R) are approximate, referring to the dominant SO interaction. We will keep 
them in the following, for convenience. 

For the bound level calculations, we extend the experimental PECS of the A and 
b states at large distances with a Cn/R^ {n = 6, 8) expansion using the coefficients 
from Ref. [50]. At large distances the functions A+(i?), A“(i?) and reach the 

SO constant of the Cs atom, i.e. = 184.68 cm“^. The constants (Abi and 
are empirically adjusted to take the relevant off-diagonal interaction into account. All 
those terms are obtained using the analytical formulas given in Refs. [45, 44, 49]. For 
completeness, we also report in Fig. lb the SOC functions which were computed in Ref. 
[51] prior to the spectroscopic analysis of Ref. [43], showing a remarkable agreement 
between the two determinations. 

The amazing quality of the spectroscopic data for the {b,A} complex allows 
calculating the energies of the O = O’*" levels with experimental precision. There are 
only a few observed levels assigned to levels with b^IIi and b^n 2 character, so that 
the prediction of the energies of unobserved levels is not as accurate. However, as it 
has been experimentally demonstrated for RbCs [52, 32], such a model indeed provides 
reliable information about these b^Hi levels. 

For the = 1 case, no full spectroscopic analysis for the b, c, and B coupled 
molecular states exists in the literature. Kim et al. computed the relevant R-dependent 
SOC functions Wbc, H4s, and Wbc with a quantum chemistry approach [51], and we 
used them in our model. The SO Hamiltonian matrix is expressed in a way similar to 
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Eq. (2) as 


V^BO + ^ = 


\h) \c) \B) 

( Vf{R) WUR) -WbB{R)\ 
WUR) Vf{R) Wbc{R) 
\-W,b{R) Wbc{R) V^iR) j 


( 3 ) 


with Vf{R) = Va{R) + fdX for a =b, c, B. The coupling functions are reported in 
Fig. 2b, showing that they all converge toward For the bound level calculations, we 
used as above the experimental b PEC [44, 49]. The c PEC is obtained from our own 
quantum chemistry calculations based on semi-empirical effective core potentials (ECP), 
following the procedure described in Refs. [40, 53], completed by the ECP parameters 
reported in Ref. [54]. The spectroscopy of the bottom of the B PEC has been achieved 
in Ref. [47], which thus accounts for SOC in an effective way. Therefore we shifted it 
in energy before its connection to our own computed PEC curve in order to ensure 
that the three-coupled-channel calculation actually delivers the correct energies for the 
measured levels. As above, these PECs are connected at large distances to an asymptotic 
expansion using coefficients from Ref. [50]. 

The treatment of the hi = 1 complex is not as accurate as the hi = O’*" one, as there is 
no spectroscopic analysis available in the literature. However, the few spectroscopically 
observed vibrational levels at the bottom of the B state provide useful information for 
predicting a good STIRAP transfer using an H = 1 intermediate level. 

In addition, TEDM functions connecting the X and a states with the A, b, B, and 
c states from our own quantum chemistry calculations are drawn in Fig. Ic and Fig. 2c. 
Note that as stated by Kim et al. in their article [51], both our own PECs and TEDMs 
are in excellent agreement with their results. It is also worthwhile to remark that 
the TEDMs are quite similar to the KRb ones, with an important difference however. 
The magnitude of df,a{R) is larger than its KRb counterpart [38] at the inner turning 
point of the a PEC: (iba(-R=9.5a.u.) ~ 0.21eao in KCs, whereas in KRb dfea(i?=9.3a.u.) 
~ 0.03eao. As discussed later, this order of magnitude difference will have important 
consequences for the experimental realization of the STIRAP approach based on D = O’*" 
states. 

Our quantum chemistry data for PECs and TEDMs and the piecewise PECs 
elaborated above are provided in the Supplemental Material for convenience. 

Finally, according to the authors, the semiempirical curves and parameters derived 
in Refs. [42, 43, 45, 44, 49] are correctly mass-invariant, so that it is possible to use them 
to model the levels of the other isotopologues ^°KCs and ^^KCs. This is partly verihed 
for the {b,A} complex by the ability of the semiempirical curves to reproduce a few 
measured levels [45] of the ^^K^^^Cs molecule. For X and a PECs, the derived curves 
yield reliable scattering lengths for elastic collisions of each isotopic combination, thanks 
to the good quality of the long range part [42, 43, 49]. 
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4. Three possible implementations of STIRAP in KCs 


In order to evaluate the relevant transition matrix elements (TMEs) involved in the 
Rabi frequencies in Eq.(l), vibrational energies and radial wave functions are computed 
with the Mapped Fourier Grid Hamiltonian (MFGH) method [55, 56] as described in 
paper I. The Hamiltonian operator governing the nuclear motion is H = T + Vbo + 1^s(P\ 
where T refers to the kinetic energy operator. It is represented as a. Nq x Nq matrix 
where N is the number of coupled channels and q is the number of grid points for the 
R coordinate determined such that the diagonalization yields eigenenergies reproducing 
the bound states energies at the experimental accuracy. The variable grid step at the 
heart of the MFGH method allows for accurately calculating bound levels with energies 
very close to the dissociation limit (and thus with a large vibrational amplitude) while 
limiting the number of grid points to about q = 590. The resulting vibrational wave 
functions |^) = | Va = 35), and \g) = | = 0), are then used in our 

calculations. 

The vibrational levels resulting from the diagonalization are labeled with an index 
Vq referring to the global numbering of the increasing eigenenergies. The corresponding 
radial wavefunctions |G; v'q) are expressed as linear combinations of the N coupled 
electronic states 

V; 4) = E l“) ■ W 

a=l 


Qv'r 


The weight Wa ^ on each electronic state is dehned by the squared radial components 


/ QVr 


f^a^(R) 


such as 


N ^ r 

0=1 0=1 “^0 


oo 2 


dR = 1 


( 5 ) 


The TMEs for the pump (resp. dump) transition involve the vibrational functions 


(resp. gfx) of a^S+ (resp. X^E+) state and the triplet part fj2t^{R) (resp. the 


, Qv'r 


singlet part ^jJaf^{R)) of the coupled wave function of the intermediate level | G; ) (eq. 

( 4 )) 


^ata 


(G, dcua I®) '^a') 

/•oo 

, Oi;; 


“Xos 


= / fjZ^{R)do^,a{R)pl^{.R)dR 

Jo 

= {X;vx\dxas 


p 7{R) dxo.XR) iS^{R)dR 


( 6 ) 


( 7 ) 


The squared matrix elements 


^ata 


and 


“Vo, 


determine the efficiency of the 


STIRAP process and are systematically calculated in the following. Note that these 
TMEs will have to be multiplied by the appropriate Holn-London factors to take in 
account the experimentally chosen polarizations of the pump and dump lasers. 
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In the next sections, graphs for TMEs will be drawn for the energy region where 
they are of comparable magnitude for the pump and dump transitions, for clarity. The 
full list of TMEs are given in the Supplementary Material attached to the present paper. 


4 . 1 . STIRAP via the A - bo spin-orbit coupled states 


We considered the lowest allowed J = 1 value for which the matrix of Eq. (2) reduces to 
& N = 3 dimension. We display in Fig. 3 for ^®KCs and in Fig. 4 for ^°KCs, the relevant 
TMEs dff^ “ (closed squares) and dxA°* (closed circles). They are extracted from the 
calculations above, involving the coupling matrix of Eq.(2). As expected, their global 
behavior is very similar for the two isotopologues, but of course the recommended levels 
for the optimal transfer are slightly different, as summarized in Table 1. 

The data points for 
eigenvectors yielded by t 


^XA 


(closed circles in Figs. 3 and 4) are associated to all 
le diagonalization, and of course related to the magnitude of 
their A component. The TMEs present strong variations associated to levels with main 
mixed A-bo character (upper zone of the data points, associated to states with strong 
A component, alternating with states with strong bo component), and to levels with 
main weight on the bi state (lower zone of the data, corresponding to states with very 


weak A component). The TMEs reported for 


u I t 

OQa 


(closed squares in Figs. 3 and 4) 


correspond to states with either a main component on bo or on A, disregarding those 
with main bi character for clarity. Thus this data is complementary to the upper part 


of the data for 


d 


XA 


Thus the optimal STIRAP region for equal pump and dump transitions, exemplified 
in the figures by the selected v'q = 202 level in ^®KCs and by Vq = 198 and Vq = 207 
in "^°KCs, is located where the upper part of the A —)■ X data crosses the a —)■ bo 
data (with closed squares). Note that the picture is qualitatively similar than the one 
obtained in paper I for KRb, except that the bi and b 2 states were not included in 
the spin-orbit coupling matrix. However, the large R-dependent TEDM around the 
inner turning point of the a PEC in KCs compared to KRb taken in similar conditions, 
namely starting from the uppermost level n = 31 of the a state (see Table 1) induces 
a much larger TME than in KRb, which makes this STIRAP scheme attractive for a 
future experimental implementation in KCs. 


4-2. STIRAP via the A - bi rotationnally coupled states 

The relevant TMEs are (open triangles in Figs. 3 and 4) and d^A^^ (fhe low set 
of closed circles in Figs. 3 and 4). The data points correspond to levels with main bi 
character and with the largest possible component on the A state (typically ~ lO”"^). 
Due to the weak rotational coupling (the constant fl in Eq.(2) amounts !^0.04 cm“^ 
around R =10 a.u.), vibrational levels of the unperturbed A and bi PECs must be 
quite close in energy to be effectively coupled. As the available KCs spectroscopic 
data is of good quality, we identified one level in ^®KCs and two levels in "‘^KCs of 








Prospects for the formation of ultracold polar ground state KCs molecules via an optical processll 


Table 1. Selected energies (in cm~^) and matrix elements (in a.u.) of the pump 
and dump transitions (with energy Epump and Edump) relevant for a STIRAP scheme 
based on an intermediate level belonging to the A - bo spin-orbit coupled states (with 
binding energy Etind relative to the 4s -1- 6p dissociation limit) resulting from the spin- 
orbit coupling between A and bo levels (with J = 1), starting from uppermost ?;a=35 
level (with J — 0). Honl-London factors are not included. The weights on the various 
components of the intermediate level are also reported. The same results are recalled 
for ^®K®^Rb [38], with as the initial level the uppermost one t>o = 31. 



v' 

Ebind 

Epump 

E^dump 

Who 

Wbi 

WA 

\dab\‘^ 

\dAx\'^ 

39KCs 

202 

-3418.866 

8128.764 

12163.843 

0.79 

0.20(-6) 

0.21 

1.21(-6) 

1.57(-6) 

40KCs 

207 

-3380.794 

8166.841 

12202.246 

0.31 

0.13^6) 

0.69 

1.40^6) 

1.40(-6) 


198 

-3500.188 

8047.448 

12082.852 

0.60 

0.17^6) 

0.40 

2.00^6) 

7.36(-6) 

39K®7Rb 

103 

-3450.1 

9287.297 

13467.396 

0.09 

- 

0.91 

9.74^9) 

2.0 (-8) 


Table 2. Selected energies (in cm“^) and matrix elements (in a.u.) of the pump 
and dump transitions (with energy Epump and Edump) relevant for a STIRAP scheme 
based on an intermediate level belonging to the A - bi rotationnally coupled states 
(with binding energy E^md relative to the 4s -|- 6p dissociation limit) resulting from 
the rotational coupling between A and bi levels (with J = 1), starting from Va=S5 
with J = 0. Honl-London factors are not included. The weights on the various 
components of the intermediate level are also reported. The results for the same 
mechanism experimentally implemented in ®^RbCs is recalled [39], starting from the 
6*^downward the dissociation limit, corresponding to Va = 42. 



v' 

E/bind 

Fj 

-^pump 

Edump 

Wfci 

Who 

WA 

Ma&P 

Max P 

39KCs 

122 

-4500.405 

7047.225 

11082.304 

0.99987 

0.00003 

0.00010 

4.42(-7) 

9.19(-5) 

40KCs 

127 

-4459.692 

7087.944 

11123.348 

0.99138 

0.00191 

0.00671 

1.37(-6) 

2.05(-3) 


162 

-3982.269 

7565.367 

11600.771 

0.99826 

0.00055 

0.00118 

1.01(-6) 

9.24(-6) 

®^RbCs 

68 

-5124.586 

6423.042 

10234.613 

0.99603 

0.00083 

0.00314 

9.49(-7) 

1.71(-4) 


main bi character with such characteristics (Table 2). Despite small TMEs and quite 
unbalanced transition matrix elements for the pump and dump transitions, we predict 
a situation which is comparable to the one modeled and already observed in RbCs, 
which characteristics are recalled in Table 2. This is actually such a circumstance which 
recently allowed for an efficient STIRAP implementation to create a dense sample of 
ultracold ®^RbCs molecules [39, 32]. This mechanism is expected to be even more 
favorable if a more deeply-bound level is chosen for ji) (as done in RbCs [39]) since 
the amplitude of its assumed pure triplet wave function around the a^S’*' inner turning 
point grows up. 


4-3. STIRAP via the B-b-c spin-orbit coupled states 

The relevant TMEs dl^Jf, dfca' ", and (fxB^ extracted from the calculations involving 
the coupling matrix of Eq.(3). They are presented in Fig. 5 for ^®KCs as a representative 
isotopologue, as the full spectroscopy of the D = 1 states is not yet available. The 
recommended levels for an optimal STIRAP implementation are displayed in Table 3. 
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Table 3. Selected energies (in cm“^) and matrix elements (in a.u.) of the pump (a —>■ c 
and a —>■ b) and dump (B —)> X) transitions (with energy Epump and Edump) relevant 
for a STIRAP scheme based on an intermediate level belonging to the B-b-c spin-orbit 
coupled states (with binding energy Edind relative to the 4s -1- 6p dissociation limit) 
resulting from the spin-orbit coupling between the b, c, and B states (with J = 1), 
starting from Uq= 35 (with J — 0). Honl-London factors are not included. The weights 
on the various components of the intermediate level are also reported. The binding 
energy of us = 0 is —1547.6 cm“^. The main vb or Vh wavefunctions involved in the 
full v' coupled vibrational wave functions are indicated. 


v' 

^bind 

Epump 

Edump 

Wb 

Wc 

Wb 

MbaP 

McaP 


39KCs 188 

{vb = 20) 

-814.4 

10733.3 

14768.4 

0.193 

0.275 

0.532 

1.32(-5) 

4.54(-5) 

3.72 (-4) 

195 
{vb = 23) 

-751.5 

10796.2 

14831.2 

0.138 

0.282 

0.579 

7.18(-7) 

6.64(-5) 

8.71(-5) 

125 
{vb = 90) 

-1543.3 

10004.3 

14039.4 

0.630 

0.357 

0.013 

3.44(-6) 

2.32(-4) 

1.08(-3) 

163 

{vb = 100) 

-1067.5 

10480.1 

14515.2 

0.528 

0.440 

0.032 

1.60(-6) 

2.83(-4) 

6.17(-4) 


As in KRb [38], due to the larger TEDM for the a —)■ c than for the a —)■ b transition, 
these levels are characterized by TMEs of comparable magnitude for the pump —)■ c 
transition and the dump B —)■ X transition. Thus we qualify this case as being induced 
by the B-c spin-orbit coupling. But in contrast with KRb, it is likely to reach quite 
high-lying B vibrational levels {vb =20, 23, while for KRb we had vb =8) to ensure an 
optimal STIRAP. These levels are located at an energy corresponding to the quantum 
chemistry part of the B PEC, while its spectroscopic determination is yielded only up 
to the energy of us = 5 in Ref. [47]. In the un < 5 energy range a couple of hi = 1 levels 
with main b character are also expected to be interesting for STIRAP and are reported 
in Table 3. For instance the v' = 125 level energy is predicted close to the location of 
Us = 0 which is well known experimentally. Due to their noticeable weight on the B 
state, such levels are most likely present in the recorded data of Ref. [47], even if not yet 
assigned. 

Finally, this B-b-c STIRAP option is based on TMEs which are larger than those 
for the A-bo and A-bi options by about two orders of magnitude only, in strong contrast 
with KRb where the difference was at least of four orders of magnitude. As anticipated 
above, this is due to the large a-b TEDM in KCs compared to the KRb one. 

5. Prospects for experimental implementation 

The present investigation of the possible pathways for the formation of ultracold KCs 
molecules in their absolute ground state is a follow-up of our previous study on KRb. 
But the current situation on the experimental side is very different for the two molecules. 

The formation of ultracold KRb molecules in their absolute ground state via a 
STIRAP scheme has been undoubtedly boosted in part by the wealth of spectroscopic 
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data available for the state in the region of its PEC minimum, with an accurate 
modeling of the perturbations by the b^II and c^S’*' states induced by spin-orbit 
interaction. Our complete analysis in paper I conhrmed the choice of the experimental 
groups for the implementation of STIRAP based on the B-b-c scheme. Indeed, the 
weak TEDM for the a^E"*" b^II pump transition does not favor the implementation 
via the other A-bo scheme relying on the spin-orbit coupled A^E’*' and b^II states. 

The situation is reversed for the KCs species which has not yet been observed in the 
ultracold regime. The spectroscopy of the coupled A^E+ and b^II states is much better 
known than the one of the B^II, b^Il and c^E"*" states. Moreover, the large i?-dependent 
TEDM around the inner turning point of the a PEC in KCs compared to KRb induces 
a much larger TMEs for the A-bo and A-bi schemes in KCs, which makes this STIRAP 
scheme attractive for a future experimental implementation. Thus our study shows 
that there are more possible options than in KRb to implement efficient STIRAP in 
KCs, namely, using either the A-bo spin-orbit coupled states, the A-bi rotationnally 
coupled states, or the B-b-c spin-orbit coupled states. Note that the A-bi case was not 
analyzed in paper I, but it has been successfully implemented in the RbCs experiment 
of Ref. [39, 32]. The magnitudes of the relevant RbCs transition matrix elements are 
recalled in Table 2, which provide a reference for the experimental feasibility of such a 
scheme in the KCs case. In this respect the STIRAP scheme based on A-bi coupled 
states is particularly attractive: it relies on the accurately known spectroscopy of the 
A-b coupled states. Moreover, despite the very weak hyperhne structure expected for 
D = O’*" levels, their coupling with closeby D = 1 levels identihed in KCs, -which are 
expected to possess a large hyperhne structure- allows for controlling the hyperhne 
level during STIRAP, as achieved for RbCs. We also predict that the STIRAP scheme 
yielding the largest TMEs for both the pump and dump transitions is the one based on 
the B-b-c spin-orbit coupled states, as observed in KRb. However, its implementation 
would require further spectroscopic investigations, which may not be the priority of the 
interested experimental groups. 

Generally, the presented study, just like the previous one on KRb, demonstrates that 
the choice of an efficient STIRAP scheme to create ultracold molecules in their absolute 
ground state level cannot be determined by invoking the Franck-Condon (EC) principle 
which only involves the spatial overlap of the vibrational wave functions. The variation 
of the relevant TEDMs along the internuclear distance plays a central role. Moreover a 
reasonable balance of the TMEs for the pump and dump transitions should be achieved, 
which generally corresponds to transitions departing from the most favorable ones 
identihed by the EC principle. In this respect, the provided Supplementary material 
should be a great help for setting up an experiment. We anticipate that the present study 
comes at the appropriate time to guide future experiments aiming at creating ultracold 
samples of KCs molecules. The hyperhne structure not taken into account in this work 
will be modeled in an upcoming study and the STIRAP schemes will reexamined in this 
framework. 
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Figure 1. (a) Potential energy curves of KCs involved in the proposed STIRAP scheme 
(schematized with thick arrows) based on the = 0“*" symmetry, and the related 
vibrational levels involved (see text). The possibility to initiate the pump transition 
from a pure singlet level is illustrated with the dashed arrows, (b) Diagonal A+ (blue 
dot-dashed line), A“ (red dashed line), and off-diagonal (solid line) spin-orbit 

coupling (SOC) matrix elements coupling the A^E+ and b^II excited molecular states 
used in the present work [44, 45]. The SOCs calculated in Ref. [51] are also displayed 
(open and closed circles), (c) Computed transition electric dipole moments (TEDMs) 
of KCs for the transition between the X and A singlet states (solid line), and between 
the a and b triplet states (dashed line), used in the present work. The same quantities 
for KRb used in Ref. [38] are displayed for comparison purpose (dotted lines). The 
TEDMs calculated in Ref. [51] are also displayed (open and closed circles). 
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Figure 2. (a) Potential energy curves of KCs involved in the proposed STIRAP 
scheme (schematized with thick arrows) based on the = 1 symmetry, and the related 
vibrational levels involved (see text), (b) Off-diagonal Wbc{R) (solid line), WtsiR) 
(dashed line), and WsdR) (dashed-dotted line) spin-orbit coupling (SOC) matrix 
elements coupling the b^II, c^E+ and B^II excited molecular states used in the present 
work [51]. (c) Computed transition electric dipole moments (TEDM) of KCs for the 
transition between the X and A singlet states (solid line), between the a and b triplet 
states (dashed lines), and between the a and c triplet states (dashed-dotted lines), 
used in the present work. The same quantities for KRb used in Ref. [38] are displayed 
for comparison purpose (dotted lines). The TEDMs calculated in Ref. [51] are also 
displayed (open and closed circles). 
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Figure 3. Squared transition matrix elements (TME) in Eq. 7 for the pump 
transitions a—i-bo (closed black squares) and a—i-bi (open blue triangles), and for the 
dump transition A—^-X (closed red circles) in the ^®KCs isopotologue, for the fl = O’*' 
case. Rotational states are: J = 1 for the intermediate levels, Jx = 0 for the hnal 
vx = 0 level, and Ja = 0 for the initial level Ua = 35. The levels Vq for which the TME 
of the pump and dump transitions are equal are indicated. 
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Figure 4. Same as Figure 3 for the "'^°KCs isopotologue. 
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Figure 5. Squared transition matrix elements (TME) in Eq. 7 for the pump 
transitions a—J-b (closed black squares) and a—i-c (open blue triangles), and for the 
dump transition B—7>X (closed red circles) in the ^®KCs isopotologue, for the 17 = 1 
case. Rotational states are: J = 1 for the intermediate levels, Jx = 0 for the final 
ujf = 0 level, and Ja = 0 for the initial level Ua = 35. The levels v'^ for which the TME 
of the pump and dump transitions are equal are indicated. 






































































